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I. INTRODUCTION 



The total nuclear-photoabsorption cross sections through the third resonance (0.2 — 
1.2 GeV photon laboratory energy) have been measured for a wide set of nuclei to study 
the behavior of baryon resonances [[lfU in nuclear matter. These measurements show an 
interesting result. In the first resonance region, the P33 resonance (I — 1, J = 3/2, / = 3/2) 
is slightly distorted. In the second and third resonance region, on the other hand, the total 
cross section is largely suppressed and the resonances such as Di 3 (Z = 2, J = 3/2, / = 1/2) 
and F 15 (7 = 3, J = 5/2, / = 1/2) disappear in the excitation function. Apparently, these 
resonances in the fundamental process, i.e., pion photoproduction off a nucleon, are strongly 
modified by nuclear medium effects. There have been several theoretical works to 
explain the disappearance of the resonances. In some phenomenological analyses |||7|], very 
large collision broadening have been assumed to fit the data. In the other work ||, however, 
it has been claimed that such significantly increasing resonance widths were hardly justified. 
The puzzle regarding the mechanism of resonance disappearance still remains unresolved. 

In the D13 and F15 resonance regions, the double pion photoproduction (7JV — > ttttN) is 
important on pion photoproduction in addition to the single pion photoproduction (7JV — > 
tiN). In those processes, the iV*(1520)(J p = 3/2",/ = 1/2) and the iV*(1680)(J p = 
5/2 + ,J = 1/2) resonances play a significant role as an intermediate state. In order to 
investigate the unknown mechanism that has caused resonances damping, one needs precise 
information about those fundamental pion photoproduction processes. 

The 7iV — > ttN reaction has been studied experimentally in the past ||. Moorhouse et al. 
Til an d Arai et al. ]lTJ analyzed the 'jN — > ttN reaction data from the first through the third 



resonance region. They made a partial-wave analysis of the processes jp — > 7i + n, ^yp — > ir°p, 
and 7n — > n~p. The imaginary parts of the amplitudes were parameterized with iT-matrices 
written as a sum of factorizable poles. The real parts of the amplitudes are calculated from 



the imaginary parts through the fixed-t dispersion relations. Furthermore, Arndt et al. [12 



made an energy-dependent partial-wave analysis on the data for the processes 7p — > n + n, 
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7P — > Tc°p, 7n — > TT~p, 7n — > 7r°n and the inverse reaction 7r _ p — > 717. From their analyses, 
D 13 and F 15 resonances have been found to be important in the 'jN — > ttN reaction in 
addition to the P33 resonance. 

The 7p — > 7t + tt~p reaction cross section has been measured as a function of the photon 
energy (0.3 — 5.8 GeV) by ABBHHM collaboration in 1968 []13|| . Recently, new improvement 
in experimental techniques make them possible to study the 7p — > tt + tt~p more accurately 
and to observe other isospin channels, i.e., 'jp — > 7r°7r + ri and 7p — > 7r°7r°p. These new data 
have been obtained at Mainz for photon energy ranging 0.45 — 0.8 GeV [14[]. Theoretically, 



the jN — > 7T7riV reactions were studied by Tejedor et al. [ TS] and Murphy et a/. ||16|| . Their 



studies show that their models can reproduce the 7p — > 7r + 7r~p reaction cross sections fairly 
well, but fail to explain the 7p — > 7r + 7r°n cross section. Tejedor et al, furthermore, calculated 
the cross sections for the neutron target, but the results are not in good agreement with the 



data In order to investigate the total nuclear photoabsorption cross section, their 

models should be improved. 

In this paper we propose a modified model for the 7iV — > ttttN reaction, taking into 
account both 7^ — > ttN and 7^ — > ttttN reactions consistently and treating the p-meson 
propagation carefully. The self-energy of the p-meson is calculated assuming the pirn form 



factor. We focus on the photon energy range measured at Mainz ||14|| . So, we include only the 
A(1232) and iV*(1520) as intermediate baryon resonance states in our model. The iV*(1520) 
will be treated carefully since it contributes importantly to the 7iV — > ttN and 7^ — > inrN 
reactions. The dominant A Kroll-Ruderman term is constructed from the finite-ranged form 
factor of the 7riVA by requiring the gauge invariance. 

This paper is organized as follows. In Sec.|| we discuss how to describe the D13 amplitude 
and the iV*(1520) propagator and furthermore how to determine strong coupling constants 
such as 7riVJV*(1520), 7rAiV*(1520) and pNN* (1520) in detail. In Sec.fTJ we discuss how to 
obtain the 7iViV*(1520) coupling constant. In Sec.[IV] we present our model of the 7iV — > 
7nrN reaction, which is based on the formalism given in Secs.|TI| and fTJ. In Sec.|V|, we show 
our predictions of the 7^ — > ttttN cross sections and the total photoabsorption cross sections 



and then we discuss the numerical results. Finally we give concluding remarks in Sec 



II. jV*(1520) RESONANCE 

The iV*(1520) resonance can decay into both the nN and tttiN channels. The branching 
fractions to ttN and -ktxN are 50~60% and 40~50% ||19|| , respectively. The tttcN decay occurs 
through three different modes, i.e., 7rA, pN and N^tt)^^^. The irA channel is in either s- 
wave or d-wave state. The branching fractions of the s-wave and d-wave decays are 5~12% 
and 10~14%, respectively. The branching fraction to pN is 15~25% and iV(7r7r)grw av e i s 
almost negligible. For simplicity, hereafter we denote iV*(1520) as iV*. 



Bhalerao et al. and Arima et al. |2T]] constructed the isobar model so as to describe 
the 7rN scattering in the D 13 channel. In this model, the N* resonance is treated as the isobar 
state. The self-energy of N* should include the contributions of both tcN and ttttN channels 



as known from the Particle Data In their models, the irirN channel was effectively 

regarded as the 7rA channel. We extend their models to apply to the pion photoproduction. 
We explicitly include three important decay modes, i.e., s- and d-wave 7rA and pN, in the 
7T7tN decay channel, since each of such modes plays a significant role in the double pion 
photoproduction process. This will be discussed later in detail. 

In the isobar model, the nN t-matrix in the D 13 channel is written as 

^fs — Mjy» — S tota i 

where yfs and M^* denote the total energy in the center-of-mass system and the bare mass 
of iV*, respectively. The total self-energy of TV*, i.e., S total , is expressed as 



J total 



SttTV + S^. A + + SpAr, (2.2) 



where S^tv, ^va^ ^va anc i ^pN are due to the coupling to the nN, s-wave 7rA, d-wave 7rA 
and pN channels, respectively. 

The vertex function for the nN — > A^* is written as 
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7t 
■kNN* 



-i(27T 



l2u>ME N (p) f, 



ttNN* 



P 



SW-Y a (pj), (2.3) 



V M ^2(m + M) \PnNN*J 
where p and p are the pion momentum and its unit- vector in the ttN center-of-mass system, 



respectively and u>^(p) = ym 2 + \p\ 2 , E N {p) = y M 2 + \p\ 2 andp = |p|. f^NN* is the itNN* 
coupling constant and p^nn* is the nNN* range parameter, and M and m denote nucleon 
and pion masses, respectively. in Eq. (|2.3| ) is defined by 



5 (2)t 



X CT 



(2) 



(2.4) 



where is the spin transition operator from 1/2 to 3/2 and cr is the ordinary Pauli spin 
matrix. 

The N* self-energy due to the coupling to the irN channel is written as 



SI 



NN* 



dp- 



p2 e -2(p/p wNN *) 2 



P 



(2.5) 



2(m + M) Jo y/s- u w (p) - E N (p) + ie \p wNN * / 
This expression is derived by using the vertex function of Eq. (|2.3| ) . 

The N* self-energy due to the coupling to the s-wave or d-wave 7rA channel is expressed 
in a similar fashion, 

d 3 p 1 



Tp S (d) 7-lS(d)t 

it AN* it AN* 



(2.6) 



(2vr)3 7u n (p) - cuM - E A (p) - E^(p, y/s) ' 
where = J (M%) 2 + \p\ 2 and is the bare mass of A. is the A self-energy due to 



the coupling to the nN channel, which expression is given in Ref. [21|. We employ the same 



7riVA vertex function and bare mass of A used by Betz and Lee |22]. The vertex functions 
for the N* — > 7rA are defined as 



7 -t 

ttAN' 



(p) = -i{2irf\ 



2uM 
\l 2(m + M) 



(2.7) 



^n(p) fd 



P 



AN* 



^AN* 



e-W*w) a (s$.Y 2 (pj), (2.8) 



\ 2(m + M)' 

where s and d denote s-wave and d-wave irA states, and /*jv*a an d f^AN* are ^ ne s-wave and 
d-wave iV* — > it A coupling constants, respectively, an* are the 7rAiV* range parameters. 
The spin transition operator from 3/2 to 3/2, sf^i m Eq. (|2.8|) is defined by 



3 , 
—m 
2 



c(2)t 
°3/2 f i 



\™>) , (2.9) 



where {jimij2m2\JM) is the corresponding Clebsch-Gordan coefficient. 
The N* self-energy due to the coupling to the pN channel is written as 

v ( n - f d3p M f pNN .fI nn . 

pN{VS> J (2n)*2u p (p)E N (p) y/8-u p (p)-E N (p) -E^fay/s)' l ' ] 

where oj p {p) = \J(rn° p ) 2 + \p\ 2 and m° p is the bare mass of the p-meson . The vertex function 



for the N* -> pN is 



Fl« N . = (2.f^*M/„,.*-)- (St . y 00 (p), (2.11) 

where e p is the p-meson polarization vector and p p nn* is the pNN* range parameter, f P NN* 
is the pNN* coupling constant. S p7r7r in Eq.( |2.1(jQ is the p-meson self-energy which is due 
to the coupling to the tttt state. The pirn vertex function is assumed to take form: 

F (m „ = 2h p (K)(e p -K), (2.12) 

where q pn7T is the pirn range parameter and f p7T7T is the p7T7r coupling constant. Using this 
vertex function, S p7r7r is written as 

SW(p, >/5) = 19 2 \ s /°° dK ^ (/lp( f ))2 (2-14) 

I2n 2 u p (p) Jo u 2 (k) ^ _ En (p) _ ^4u;2( K ) +p 2 + ie 



where f K = 2uj. k (k) / y 4uj 2 (k) + p 2 . The above expression for the p-meson self-energy is 
obtained by extending the self-energy in the rest frame of p to that in the moving frame. It 
should be noted that the imaginary part of Eq.( ^TJ ) has a right form of the half width for 
the p-meson with the momentum p, 

= 9 / n f , {h p {n e )f , (2.15) 

2 24:7TUJ p {p)UJ n (K e ) 

where n e satisfies yfs — E^(p) — ^Auj 2 (K e ) + p 2 = 0. When m° p is replaced by the on-shell 
p-meson mass and q pn7T is taken to be infinite, Eq. fl2.15p is reduced to the width of Ref. |15 



In order to calculate the self-energy of the p-meson, one needs to know three independent 
parameters, i.e., f p7rn , q pn7T and m° p . These parameters can not be determined uniquely by 
using the mass and width of the p-meson. Therefore, we treat the parameter q p7T7T as a free 
parameter and vary it to reproduce the 7p — ► 7T + 7i°n data. If q p7Tn is fixed, m° p and f p7Tn are 
determined by the following condition, 



m e p xp - 



(m° + S p _(0,m e p x P + M)) » z^ [MeV], (2.16) 



where mf p « 770 MeV. 

The self-energies and bare masses of A and p in the 7V*-propagator are obtained through 
the data such as the ttN P33 scattering, the resonance energies and their widths in a phe- 
nomenological way as mentioned above. So there are nine parameters which have to be 
determined: the coupling constants (f n NN*, fnAN*> fnAN*> fpNN* ) > the range parameters 
(p-irNN*, PnAN*i PnAN* > VpNN*) and the bare mass (Mjy>). We determine the parameters by 
fitting them to branching ratios, the N* resonance energy, its width and the energy depen- 
dence of the ttN D13 scattering amplitude. In our model, we use 1520 MeV as the resonance 
energy and 120 MeV as the width. We take a fraction of 58% for the decay into irN, 10% 
into s-wave 7rA, 10% into d-wave 7rA and 22% decay into the pN channel, respectively. 

The parameters obtained are given in Table |l[ The parameter-set (I), (II) and (III) have 
been obtained by using Eq. (|2.12| ) with q pnw = 00, 100 MeV/c and 200 MeV/c, respectively. 
In the parameter-set (I), the self-energy and the bare mass of p in Eq. fl2.10p have been 
assumed to be the width and the on-shell mass, respectively, which corresponds to the 



treatment of Ref. ||15|| . One finds that the ttN D 13 partial- wave amplitude can be equivalently 
reproduced by any parameter-set in Table |. It should be noted that the sign of coupling 
constants is not determined from the data, since coupling constants appear as their squared 
form in the self-energy. Which parameter-set and which sign are appropriate will be discussed 
in Sec.|V|. 
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III. jNN* COUPLING 



In this section, we show how to determine the <yNN* coupling constants. Obviously, the 
ryNN* vertex has two independent helicity couplings. For the proton target, the helicity 1/2 
amplitude is small enough compared with the helicity 3/2 amplitude. Hence, helicity 1/2 
amplitude could be neglected |PH|-[T^]. For the neutron target, on the other hand, one can 
not use this approximation, since the helicity 1/2 amplitude is non-negligible. 

The resonant amplitude in the isobar model has the form 



T$ - F nNN *— — F NN ,, (3.1) 

y/S - M N * - Ltotal 



where F^nn* is the vertex function for the 7 A" — > N* transition. For the helicity 1/2 
transition, F^ NN , is written as 

F%* = -Wi/2 (S j -k)(a-kx e), (3.2) 

where gi/2 and £ are the helicity 1/2 coupling constant and photon polarization vector, 
respectively, and k denotes the unit vector of initial photon momentum. In our approxima- 
tion, this helicity coupling constant for the proton is set to zero. The helicity 3/2 transition 
F^ NN * operator is written as 

F%%* = g 3/ 2 {0 f • e) + \ {& ■ k) (cr ■ k x s)} , (3.3) 

where c/3/2 is the helicity 3/2 coupling constant. 

The relevant multipole amplitudes can not be described by only the resonant form of 
Eq. ( |3.1|) since there is a non-negligible background process. Actually, the full D 13 amplitude 
should be expressed as the sum of two terms, i.e., the background and the N* resonant 
terms (Fig© j^H, 

T(D 13 ) = T B + f N *, (3.4) 

m _ F I pl/2,(3/2)f / o K\ 

y/S - M N * - Ltotal 
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where Tb is the non-resonant multipole amplitude which is obtained from the partial-wave 
decomposition of the Born term in the 7 iV — > ttN reaction. The jNN* vertex function is 
rewritten as 

p l/2(3/2)f 

£,1/2(3/2)1 ~ / r-\ r 'yNN* (o a \ 

F jNN* =91/2(3/2) Ws)— , (3.6) 

#1/2(3/2) 

where #1/2(3/2) represents the effective helicity 1/2 (3/2) coupling constant, which includes 
the contribution of the N* production through the irN, irA and pN intermediate states 
(Fig.Q (c) to (e)) in addition to the direct production (Fig. |^ (b)). This coupling constant 
is complex and energy dependent. 

The effective jNN* coupling constant #1/2(3/2) is phenomenologically determined by 
through a fit to the experimental multipole amplitudes of Refs. [[TD|^T^] , instead of cal- 
culating the diagrams of Fig.^| (c),(d) and (e) in a microscopic way. Here we use the Born 



term with the same coupling constants and cutoff employed by Nozawa et al. p4| . If the 
multipole amplitude of Ref. [12|] is used, for example, the helicity 3/2 coupling constant for 
the proton target becomes #3/2 = 0.1621 + z0.0522 at 750 MeV photon energy. 



IV. MODEL OF THE jN — ► vrvriV REACTION 

The total cross section for the 'yN — > -k 1 -k 2 N reaction is given by 

1 Ml f d 3 p f d 3 q 1 d 3 q 2 M 1 1 
a ~2kE~V s l (2vr) 3 (2vr) 3 (2vr) 3 T f 2u n { gi ) 2^(g 2 ) 

x (2tt) 4 5^( Pl + k - p f - q x - q^ W '\ \ (1/2, v \T\ 1/2, v')\ 2 , (4.1) 

where pi = (£/j,pJ, pf = (Ef,Pf) and q a = (u n ,q a )(a = 1,2) are the initial nucleon, the 
final nucleon and the final pion (tt 1 ' 2 ) 4-momenta in the center-of-mass system, respectively, 
and v s is the relative velocity of the initial nucleon and the photon. The absolute square 
of the invariant matrix element T for the 7iV — > ttttN reaction is summed over the final 
nucleon spin states (z/) and averaged over the initial nucleon spin states (y). 

We describe how the matrix element T in Eq.( |4.1| ) is derived within our approach. We 
assume that the 7^ — > ttttN reaction is dominated by the processes of the 7 iV — > irA — > imN 
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and the 7 AT — > N* — > ttttN. In this assumption, there are four important processes shown 
in Figs. 0(a), (b),(c) and (d). We neglect other possible diagrams involving A, which are 
obtained from the requirement of the gauge invariance, since these contributions has been 
shown to be small | fL5|j . 

The diagram (a) in Fig.|3] contains the ^NttA contact term, i.e., the A Kroll-Ruderman 
term. This ■yNirA contact term operator -F^kr is obtained from the strong uNA vertex 

we 



function by requiring the gauge invariance. Instead of using the effective Lagrangian [T5| 
start from the vertex function with a form factor. The N — > nA transition operator F^ NA is 
assumed to be the same form with the A — > ttN vertex function which is phenomenologically 



given in Ref. [Q. Since the range parameter (see Appendix) may not be necessarily 
the same, we treat it as a free parameter and vary it to fit the 7£> — > -K + 7r~p cross section 
(see Table |). Its matrix element in coordinate space in non-relativistic limit is given by 
(suppressing the isospin factor) 



7T 



A 



r irAN 



N 



- V p )H(p)* N (r N ), (4.2) 



where r^- and p are the 7rA center-of-mass and relative coordinates, respectively, and 
MjrA = + m. The ttNA form factor H is given in Appendix. The A and pion wave 
functions in Eq.( [4.2[ ) may be expanded around r^r in power series of the relative coordinate 
p. If all gradients with respect to operating on the A and pion wave functions are re- 
placed by V ' n — ie^A(rM) or Vat — ze 7r A(rAr), where e& and e n are the A and pion electric 
charges, the resulting vertex function is invariant under the gauge transformation p5j . The 



electromagnetic interaction for the 7^ — > ttA process is then derived by the expansion to 
order A. Simplifying this interaction further according to the prescription described in Ref. 



2q| , we obtain 

fIkk = -1 {Gi{& ■ e) + G 2 (S^ ■ q a )(q a ■ e)} . (4.3) 

This expression, i.e., the minimal interaction current, is employed in our model. Here, masses 
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of the pion and the A in Eg. ( |4.2|) are replaced by their energies as the relativistic general- 
ization and the jNirA form factors Gi, G 2 are defined in Appendix. The first term in the 
left-hand side of Eq. Q4.3D corresponds to the ordinary contact term (the A Kroll-Ruderman 
term). The second term, on the other hand, appears due to the presence of the 7riVA form 
factor. It should be noted that the above minimal interaction current has a contribution 
to the 7p7r°A + vertex because of the charged A, although there is no contribution within 
the framework of the effective Lagrangian. As shown later, this interaction has a small but 
non-negligible effect on the ^p — > 7r°7r°p cross section. 

In our model, the 7./WA pion-pole term F^ PP included in Fig.^(b) will be derived from 
the time-ordered perturbation theory. The vertex functions for the A — > nN and iV — > 71 A 
transitions are of the same form as that in Eq.( f4.2|) . But the range parameter Q A of the 
latter is taken to be the same as that of Eq. (|4.3|) . The operator F APP is then written as 

APP 2^{\q a -k\)\ D A (q a , E N (k) - u„(\q a - k\)) [ l) 

H{ "' 2) [&. K2 ))( qa .e), (4.4) 



k - u n (q a ) - tUn(\q a ~ k\) 
where 

E A {q q )k 

E A {qa) + Undid - k \) 

- 1, E N{k)q a . s 

" 2 " E N (k)+^(\q a -k\Y (4 ' 6) 
D A (p,E) = E- E A (p) - rtl N) (p,E). (4.7) 

Here, Ki and k 2 are the 7rA and ttN relative momenta in the intermediate state and 
E ^ ^ {p, E) is the self-energy of A with the momentum p. The charge-dependent factor 
g p is given in Appendix. 

Using Eqs.( [4.3|) and ( f4.4|) , we can write the invariant matrix element for the diagrams of 
FigfJ(a) and (b) in the following. 

rp _ ^jVA^AKR.(PP) , 

Jakr,(pp)-— — — — -^y — , (4.8) 

y/s -u n (q a ) - E A (q a ) -T, A '{q a ,V s > 
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where Takr and Tapp represent the A Kroll-Ruderman term and A pion-pole term, respec- 



tively. Multiplying Eq.([4.8|) by appropriate isospin factors, we get the T-matrix elements 
for various reactions. 

For the 7iV — > N* — > ttttN process, there are two possible processes accompanied with 
either the N* — > ttA or N* — > pN decay as shown in Figs.|3](c) and (d). The N* resonance 
can decay into both s-wave and d-wave 7rA states. In order to construct the T-matrix 
involving the iV*, we use the strong and electromagnetic vertex functions obtained in Secs.|Tl 
and 0. Using Eqs.Q, Q, (|2~Tl| ), and (gj), the T-matrix elements of Fig. |(c) 

and (d) are written as 

TP 77 s (0t pt 

^s(d)-wave _ r nNA-F ^AN**^ ~{NN* gx 

N ' (7TA) ~~ Ws - u>M - E A (q a ) - Z% N) (q a , y/5)) - M N * - S total ) ' 

N * ipN) 2u p (q p ) (Vi - ojMp) - E N {q p ) - S p _(g p , - M N * - E total ) ' 1 ' J 

respectively. Here, q p = \q 1 + q 2 \. 

Therefore, the invariant matrix element T in Eq.(|4.1|) is expressed as follows: 



rp rp | rp , rps — wave , rpd — Wave I rp / /I 

-L - -L AKR + J APP + -tAT*(7rA) + 1 V*(ttA) + 1 N*{pN)- l 4 - 1J -J 



V. NUMERICAL RESULTS AND DISCUSSION 

In this section, we present our calculations of the total cross sections for the 7 A" — > imN 
reaction, which are shown in Figs.|] to || We calculated them with the parameters obtained 
in previous sections, but the sign of some strong coupling constants and the range parameter 
Qa{N — > 7rA) will be determined so as to fit them to the ^p — > tt + tt^p data. Then, we will 
discuss which parameters of the pirn form factor are favored by the jp — > 7T + 7i°n data. In 
our numerical calculations, the Monte Carlo integration package BASES25 |?B] is used. 

At first, we show the results of 7p — > 7T + 7T~p and 771 — > 7T + 7i^n reaction cross sections 
(solid lines) in Figs.^-(i) and (ii) which are calculated with the parameter-set (I) in Table 
|. As can be seen from Fig.|]-(i), the A Kroll-Ruderman term Takr and A pion-pole term 
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Tapp terms (dashed line) dominate on the ^p — > tt + tt~p reaction. We observe that the 
N* contribution (dash-dotted line) alone is small but the interference between the N* term 
^n~*JwA) an d T&kr is important. This strong interference occurs due to the fact that T^~J^ 
has the same structure as Tkr. Because of this, the N* excitation is regarded as an important 
ingredient in the 7J9 — > itttN reaction. 

As mentioned in Sec.|J, there is an ambiguity about the sign of strong coupling constants, 
fnAN* an d fpNN* • We adopt a positive sign for f*& N * which gives rise to a constructive 
interference between T^~*J^ and Takr- The peak position of the jp — > n + n~p cross section 
can be reproduced with this choice as shown in Fig.^-(i). We find that the curve with 
fwAN* > 0; fnAN* < 0; an d fpNN* < agrees well with the jp — > 7t + ii~p data |13| , |T^,[l7| , [18 



The range parameter Qa{N — > ttA) is taken to be 420 MeV/c. Hereafter, we will use the 
same sign for the coupling constants. Furthermore, the — > 7i + 7i^n reaction cross section 
(Fig.|-(ii)) is calculated by using the same parameter-set (I) except for the ^NN* coupling. 
Our calculation with this parameter-set is also in good agreement with the jn — ► -K+n^n 
data P7|jl5| . This is different from the result by Tejedor et al. [13]. Their model could 



not reproduce the — > 7i + 7i^n data in spite of good agreement with the 'jp — > tt + tt~p 
data. This difference may be attributed mainly to the theoretical treatment of the A Kroll- 
Ruderman term. In any case, new experiments for the neutron target would be welcome in 
order to check the validity of our model. 

Secondly, we show the results of the 7p — > 7i + 7i°n and — > 7t~tt°p reactions in Figs.^-(i) 
and (ii) (thin-solid lines), which are calculated with the parameter-set (I). In both reactions, 



we find a large discrepancy with the data ||14||i7||18f . Our model underestimates cross sections 



about a factor of two compared with the data. In these reactions, as shown in Fig.[5|-(i) and 
(ii), Takr and Tapp contribution (short-dashed line) is very small. We also find that the 
contribution of the N* terms (dash-dotted line) is almost the same as that in the jp — ► -k + -k^p 
reaction but there is no characteristic energy- dependence due to the interference which 
is clearly observed in the 7^ — > -k + -k^p cross section. This can be understood by the 
isospin factor. The isospin ratio of the Takr(pp) term of the 7p — > 7i + 7T°n reaction to the 
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7£> — > vr + 7r p reaction is 

T A KR(PP)(7T + A -> 7T + 7T r2) _ _ 



7AKR(PP) (tT A++ -> 7T 7T+p) 3 

where only the intermediate and final states are written in the invariant matrix element 
T but the initial state 'jN is omitted. For instance, T(tt~A ++ — > 7r + ir~p) means T^p — > 
7r~A ++ — ► 7T + 7r _ p). Eq.( |5.1|) indicates that the cross section for the jp — ► 7r + 7r°n reaction is 
2/9 times smaller than 7£> — > 7r + 7r~p reaction. For the A Kroll-Ruderman term, this relation 
holds exactly if only the dominant term is considered, but this feature is not changed even if 
other small terms are included. The isospin ratio of the T^f^~^ ve term of the 7p — > ir + TT°n 
to 7p — > 7t + tt~p reaction is 



Eqs. (|5.1| ) and (|5.2Q show that a relative sign between T K R(pp) and T^lJ^ in the 7p 



(5.2) 



7i + n n 



reaction is different from that in the 7p — ► -k + tt~p reaction, and Eq.( |5.2|) indicates that 
^n^OkA) m ^ ne TP ~~ * 7T + 7r°?T- reaction is smaller than that in 7£> — > 7r + 7r _ p reaction. Therefore, 
the interference between Tkr and T^lJ^ in the 7£> — > 7i + 7i°n reaction is different from that 
in the 'jp — > 7r + 7r~p reaction. In addition, there is another important feature regarding the 
p-production amplitude (diagram (d) in Fig||). The isospin ratio of the T N *( pN j term is 

The Eq. (|5.3| ) shows that the Tn*( p n) term in the 7£> — > 7r + 7r°n reaction is larger than that 
in the 7p — > tt + tt~p reaction. Hence, the p-meson production term is important in the 
7P — > 7r + 7r°n reaction. The large shift of the peak compared with the calculation of the 
7P — > n+n^p cross section may be due to the large T N */p N \ term. The same arguments 
remain true for the — > 7r + 7r~n and 'jn — > 7r~7T°p reactions. 

As can be seen from Figs.^-(i),(ii), Figs.^J-(i) and (ii) (thin-solid lines), one finds that our 
model is successful in the 7iV — > tt + tt~N reaction but fails to reproduce the experimental 
results of the 7p — > 7r + 7r°n and — > 7T~ir°p reactions. This difficulty concerning the ir° 
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production has been already pointed out in other studies |T5|JT6[] . In order to improve our 
model and resolve this problem, we examine the two possible processes shown in Figs.|3] (e) 
and (f), which may contribute more effectively to the jp — > 7T + 'K n than the jp — > 7r + 7r~p 
reaction. The diagram (e) in Fig.^ is the final state rescattering process. But we found from 
our rough estimate that the cross section for this rescattering diagram (e) is very small. 

Next, let us discuss the effect of the diagram (f) in Fig|| in detail. This diagram con- 
tributes only to the 7p — > 7i + 7T°n and — > 7T~7r°p reactions, since it contains ^yNpN contact 
interaction (p-Kroll-Ruderman term). The invariant matrix element for this diagram is writ- 
ten as 

TP /?t 

pKR 2u p (q p )(y/s - uj p (q p ) - E N (q p ) - E p7T7T (q p , y/s))' 
Here F^ NpN represents the •yNpN contact term (suppressing the isospin factor) 



F^npN = iefc\j ENi 2M M ^M Sp ' ( °" X 6) ' (5 ' 5) 



where G T = —17.6 is the tensor coupling to the pNN channel and the factor f c is taken to 



be ^Jm^rn'p/ (m c p xp ) 2 . m' p = m° p + S p7r7r (0, M) is the mass of p carrying zero momentum and 
zero energy. The factor f c arises from the fact that the intermediate p-meson mass in our 
model is different from the on-shell mass. The effect of this factor is however small so that 
it does not affect our final results significantly. Thus, our new invariant matrix element T 
matrix becomes 

Trp I rp I rps— wave i rpd— wave , rp , rp a\ 

— -L AKR + J- APP + i- n*(ttA) + 1 7V*OA) + - Z A r *(pA r ) + 1 pKR- [o-O) 

We neglect the other possible diagram obtained by requiring the gauge invariance, i.e., the 
p-meson pole diagram, since its effect was found to be negligible. 

In order to see the effect of the diagram (f) in Fig.^, we calculated the jp — ► 7T + 7i°n 
cross section by using the invariant matrix element T of Eq. (|5.6|) , which is shown as dotted 
line in Fig.^|-(i). In this calculation, we used the parameter-set (I). As far as T p kr term is 
concerned, in this case, the theoretical treatment in our model is essentially the same with 
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the model by Murphy et at [|T^]. From the comparison between the calculation with T p kr 
(dotted line) and without T pK R (thin-solid line), the contribution of the diagram (f) in Fig.|3| 
is found to be very small, which is consistent with the result of Ref. [[UJ . 

To improve our model, we will treat the p-meson propagator in a dynamical way since the 
p-meson involved in the double pion photoproduction is not on-shell below the energy range 
800 MeV. We allow the range parameter of the p7T7r form factor to be finite as described 
in Sec. ||. We calculate the cross section with the parameter-set (II) in Table H where the 
range parameter of the pirn form factor g p7r7r is 100 MeV/c and Qa{N — > irA) is 400 MeV/c. 
For the 7p — > tt + tt~p cross section, the result with the parameter-set (II) is almost the same 
as that with the parameter-set (I). As can be seen from Fig.|5|-(i), however, the significant 
enhancement occurs in the cross section of the 7p — > 7c + 7c°n reaction (bold-solid line). Our 
improved model fairly well reproduces the data of the 7p — > 7i + 7i°n reaction, except for 
the energy region above 750 MeV. This result is quite different from those by Tejedor et 
at [plj and Murphy et at ||16||. Furthermore, we have calculated the — ► 7r _ 7r°p reaction 



cross section by using the same parameter-set. The calculation is shown in Fig .|§-(ii) (bold- 
solid line). We find that our model with this parameter-set is also able to reproduce the 
experimental data for the neutron target |]TTL[ 



In order to find the reason why such significant enhancement has occurred, we plot a 
27r-spectral function, i.e., ImS p7r7r /|A/i — m p — S p7r7r | 2 which is proportional to the integrated 
cross section for the diagram (f). The curves with the range parameter q pn7T = 100 (solid 
line), 200 (dashed line), 300 (dash-dotted line) MeV/c, respectively, are plotted in Fig.[7|. 
The peak is clearly seen at the p-meson resonance energy. Furthermore, below the 400 MeV, 
there is a small bump in the curve with q pn7T = 100 MeV/c. This non-negligible effect may 
cause a significant enhancement to the 7p — > 7T + 7T°n reaction cross section. One can see the 
similar enhancement in the 27r-spectral function associated with the NN — ► 2n reaction in 
Ref. ||27| . This low-energy behavior comes from the background contribution in the isospin 
1 = 1 channel. Our model for the 2% scattering includes only the pirn coupling, but not the 
background interaction. However, the small range parameter in our model might simulate 
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the background interaction effectively. 

We show the results of double neutral pion photoproduction processes in Figs.|6|-(i) and 
(ii). As can be seen from Fig.^-(i), the peak position of the calculated cross section of the 
7P — > 7r°7r°p is in good agreement with the data On the other hand, the magnitude of the 
cross section is underestimated about a factor of two. Murphy et al. suggested |TBj that the 
7P — > TT°n°p reaction cross section was enhanced by the presence of the Pn(1440) resonance 
which decays into the cr-meson. However, it seems that other important mechanisms are 
still missing. We leave this problem as a next step since we are interested in the total 
photoabsorption at present. We also show the 772 — > 7r°7r°n reaction cross section in Fig.^- 
(ii). Unfortunately, there are no experimental data to compare with our calculation. 

Finally, we show the total photoabsorption cross sections (solid lines) for 7p and 772 
reactions in Figs.§-(i) and (ii), respectively. We assume that the total photoabsorption off 
a nucleon is dominated by the 7iV — > 71N and 7iV — > ttttN reactions in the energy region 
where we discuss. The 7iV — > ttN cross sections (dashed lines) are calculated by using the 
amplitude given in Ref. [T2|] except for the D 13 amplitude. The D 13 amplitude is treated in 
the same way described in Secs.y and fVj. On the other hand, the 7iV — > itttN reaction 
cross sections (dash-dotted lines) are calculated by using our model with the parameter- 
set (II). We found our model reproduces the experimental data both ■yp [U^U an d 7^ |29| 
reactions over the wide range of the energy. We should mention that the double neutral 
pion photoproduction scarcely contribute to the total photoabsorption cross section. 



VI. CONCLUSION 

We have constructed a simple model for the 7iV — > tittN reaction. It is assumed in this 
model that the processes of 7^ — > 7rA(1232), 7^ — > iV*(1520) and 7iV — * pN are dominant 
in the double pion production. We treat the resonances such as A, iV* and p-meson in a 
dynamical way. 

The 7iWA contact operator is derived from the strong 7riVA vertex function by requiring 
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the gauge invariance, instead of the effective Lagrangian. The range parameter of the N — > 
7rA form factor in this operator is determined so as to reproduce the 7p — > -k + -k~p data. For 
the A resonance, the dynamical model by Betz and Lee [^2| was used. We have carefully 
discussed the dynamical model for the N* resonance because this resonance plays very 
important roles in the energy range where we concern. The strong vertex functions of 
nNN*, ttAN* and pNN* are obtained from the nN scattering amplitudes, decay widths 
and the resonance energy of N*. The sign of each strong coupling constant is not determined 
from only the above data. However, as the sign affects the energy dependence of the double 
pion photoproduction cross section, it could be fixed to reproduce the 7p — > 7r + 7r~p data. 
The electromagnetic couplings of ^yNN* are determined through a fit to the 7iV — > tcN D 13 
helicity 1/2 and 3/2 partial- wave amplitude. 

Our model with the above parameters can simultaneously reproduce total cross sections 
of both the ^p — > 7T + 7T~p and — > 7i + 7i^n reactions. We found that the A Kroll-Ruderman 
term and the A pion-pole term had a dominant contribution to these reactions and the 
interference between the A Kroll-Ruderman term and the N* term was very important to 
explain the energy- dependence of the cross sections. 

For the 7£> — > 7r + 7r°n and 'yn — > 7t~tt°p reactions, on the other hand, the contributions 
of the A Kroll-Ruderman term and the A pion-pole term are smaller compared with the 
7iV — > 7T + 7T~N reaction. In order to reproduce the 'jp — > 7r + 7r°n data, we found that one 
should treat the p-meson in a dynamical model where the pirn vertex has a finite-ranged form 
factor. In fact, it has turned out in our calculation with the appropriate range parameter that 
the syNpN contact term contributes to these tt° productions significantly. As a result, our 
improved model could simultaneously reproduce both the 7p — > 7i + 7i°n and 'yn — > tx~-k°p 
cross sections, except for the energy region above 750 MeV. The disagreement at higher 
energy may be due to the fact that other higher resonances are not taken into account in 
our model. 

The total cross section of the ^p — > TT ir°p reaction is underestimated about a factor of 
two. Still some important mechanisms are missing in our model. We leave this problem as 
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a next step. 

Finally, we have calculated the total cross sections of the photoabsorption off a nucleon. 
We found that our model is able to reproduce the experimental results of the proton target 
as well as the neutron target. The defect of our model regarding the double neutral pion 
production gives little influence on the total photoabsorption cross section, since the mag- 
nitude of its cross section is very small. As far as the total cross section is concerned, we 
consider that our model has a predictable power for nuclear processes. Based on our model, 
we are investigating the mechanism which has caused the resonance damping in the nuclear 



photoabsorption ||30|| . 
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APPENDIX A: THE vriVA AND 7 iV7rA FORM FACTORS 

The irNA form factor H is defined as 



H(q) = ^ ^MEMg^jq) ^ (A1) 



where q is the 3-momentum in the tiN center of mass system and g^NA is given by [22 



q ( Ql V 
^/2(m + M) m \Qa + Q J 

where is the coupling constant and Qa is the range parameter. The 'jNnA form factor 

G\ and G 2 in Eq.( [4.3|) and the factor g p in Eq.(O) are given as follows. 

1) 7p — > 7r~A ++ reaction 

G x = e {Z 2 H(\q - Z 2 k\) + 2Z 3 H(\q + Z 3 k\)} , (A3) 

G 2 = e{-~h n + 2~h A ), (A4) 

g P = -2e. (A5) 
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2) 7£> — > 7T + A° reaction 



G 1 = --^Z 2 ^(|q-Z 2 fc|), (A6) 

G 2 = (A7) 
2 

<? P = ~jf- (A8) 



3) 7p — > 7r°A + reaction 



G 1 = --^Z 3 ff(|g + Z s fc|), (A9) 

G 2 = eh A , (A10) 

<?p = (All) 

4) — > 7r + A~ reaction 

Gi = e {-Z 2 if (|g - Z 2 fe|) - Z 3 H(\q + Z 3 k\)} , (A12) 

G 2 = e(h 7T -h A ), (A13) 

£ P = 2e. (A14) 

5) 7n — > 7r~A + reaction 

Gi = ^= {Z 2 H(\q - Z 2 k\) + Z 3 H(\q + Z 3 k\)} , (A15) 

G 2 = -^(/ia-M, (A16) 

9p = ~^=e. (A17) 

where e is electromagnetic charge and including isospin Clebsch-Gordan coefficients of the 
jNttA vertex and q is the 3-momentum of the out-going pion and 

Z2 = T-i E t iq \^ (A18) 

Zs = " (?) , r , (A19) 

K = ^T^-^ , (A20) 
k - (q- Z 2 k/2) 

~ H(\q + Z 3 k\)-H(q) 

= k.(q + Z 3 k/2) ■ (A21) 
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TABLES 

TABLE I. The parameters used in our model. Only the absolute value of coupling constants 
is given. The sign is discussed in the text. 





parameter-set (I) 


parameter-set (II) 


parameter-set (III) 


M N * (MeV) 


1597 


1554 


1566 


fnNN* 


1.09 


1.13 


1.13 


P-kNN* (MeV/c) 


450 


400 


400 


fnAN* 


0.992 


0.992 


0.992 


PnAN* (MeV/c) 


200 


200 


200 


fwAN* 


0.984 


1.00 


1.00 


p d wAN * (MeV/c) 


200 


300 


300 


fpNN* 


1.56 


0.928 


0.583 


PpNN* (McV/c) 


200 


200 


300 


fpmr 


6.14 


82.0 


25.6 


q pnn (MeV/c) 


oo 


100 


200 


Q A (N -» ttA) (MeV/c) 


420 


400 


400 


Qa(A -» ?riV) (MeV/c) 


358 


358 


358 
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FIGURES 

FIG. 1. The 7iV — > irN reaction in the D13 channel including background and iV* production. 
B.G. corresponds to the background term, (a) The full-Di3 amplitude, (b) The Born term for the 
jN — > irN reaction, (c) The N* resonant term through the effective jNN* vertex. 

FIG. 2. The effective -yNN* vertex, (a) The effective -yNN* vertex, (b) The bare 'jNN* 
vertex, (c) The vertex correction due to background irN production, (d) The vertex correction 
due to background it A production, (e) The vertex correction due to background pN production. 

FIG. 3. The diagrams for the —* nirN reaction, (a) The A Kroll-Ruderman term, (b) 
The A pion-pole term, (c) The N* — > ttA contribution, (d) The N* — > pN contribution, (e) The 
final state rescattering contribution, (f) The p-meson Kroll-Ruderman term. 

FIG. 4. Total cross sections for (i) the jp — ► Tr + ir~p and (ii) the — > 7r + 7r~n reactions. 
Solid line corresponds to the total cross section calculated with the parameter-set (I), dashed line 
to the contributions of the A Kroll-Ruderman and A pion-pole terms (diagrams (a) and (b) in 
Fig|3|) and dash-dotted line to the contributions of the N* terms (diagrams (c) and (d) in Fig|3|). 
Experimental data are taken from Refs. p^,14,17,I8[. 



FIG. 5. Total cross sections for (i) the 7^ — > 7r + 7r°n and (ii) the 772 — > tt tt°p reactions. 



Thin-solid line corresponds to the total cross section calculated by using Eq.( 4.11 ) with the pa- 



rameter-set (I) and dotted line to the total cross section calculated by using Eq.(|5.6|) with the 



parameter-set (I). Bold-solid line corresponds to the total cross section calculated by using Eq.(5.6) 
with the parameter-set (II) , short-dashed line to the contributions of the A Kroll-Ruderman and 
A pion-pole terms (diagrams (a) and (b) in Fig.^), dash-dotted line to the contributions of the 
N* terms (diagrams (c) and (d) in Fig||) and long-dashed line to the contribution of the p-meson 
Kroll-Ruderman term (diagram (f) in Fig]3|) including the finite-ranged form factor of pirn. Ex- 
perimental data are taken from Refs. p^Jl7|,p!8|. 
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FIG. 6. Total cross sections for (i) the 7p — > -K ir p and (ii) the — ► 7r°7r°n reactions. Solid 
line corresponds to the total cross section calculated by using Eq. (|5.6|) with the parameter-set 
(II), dashed line to the contributions of the A Kroll-Ruderman term (diagram (a) in Fig.||) and 
dash-dotted line to the contribution of the N* term (diagram (c) in Fig.||). Experimental data is 
taken from Ref. 14]. 



FIG. 7. The 27r-spectral function as a function of the center-of-mass energy. Solid line corre- 
sponds to q pnw = 100 MeV/c, dashed line to q pn7T = 200 MeV/c and dash-dotted line to q pn7T = 300 
MeV/c. 

FIG. 8. The total photoabsorption cross section of (i) the proton target and (ii) the neutron 
target. Solid line corresponds to the summed cross section of 7iV — ► irN and 7iV — > 7T7T./V which are 
calculated in our model. Dashed line corresponds to the contribution of 7^V — > irN, dash-dotted 
line to the contribution of 7./V — * ttttN. Experimental data are taken from Ref. jj,28,2Jj. 
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